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Abstract 

The inverse of the Vandermonde and confluent Vandermonde matrices are pre- 
sented. In the case of the Vandermonde matrix, we present a decomposition in three 
factors, one of them a diagonal matrix. The evaluation of such inverse matrices is a key 
point to find functions of a matrix, namely exponential functions (evolution operators) 
and logarithmic functions (entropies) in quantum mechanical topics. 
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1 Introduction 

Although Vandermonde systems arise in many approximation and interpolation problems 
[TJ [5] , they also appear when we need to solve systems of differential equations [3] , such as 
when we interact a multi-level atom with a classical or quantum field, a trapped ion with 
a laser field (see for instance [1]), etc. When an atom interacts with a quantized field they 
get entangled |5J. This produces that, by analyzing the density matrix of the atom or the 
density matrix of the quantized field, we can determine when they disentangle. To do so 
we need either to calculate the entropy of the sub-systems and this requires evaluation of 
functions of (density) matrices. 

In the particular case of the entropy, we need to calculate logarithmic functions of the sub- 
system's density matrices [S]- Vandermonde matrices, and in particular, their inverse, are 
helpful to determine such functions. A more common function is the exponential function 
of a matrix, as a Hamiltonian may be written usually in matrix form, and therefore the 
solution of Scrodinger equations involve the use of evolution operators, i.e. exponentials of 
Hamiltonians (see for instance [B], also for the case when superoperators are considered, and 
[7] for time dependent Hamiltonians). The key point for the evalution of such functions is 
to find the inverse of a Vandermonde matrix or of the confluent Vandermonde matrix (in 
case there are repeated eigenvalues). The purpose of the present paper is precisely this. 
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2 Vandermonde matrices 



A matrix N x N of the form 
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Vi,i = A}" 1 



i = 1,2, 3, ...,N; j = 1,2, 3, ...,7V 



(2) 



is said to be a Vandermonde matrix [8j [T] . 

The determinant of the Vandermonde matrix can be expressed as 

det(tO= II (Ai-Ai). 

l<i<j<JV 

Therefore, if the numbers Ai, A2, \n are distinct, V is a nonsingular matrix [8]. 
When two or more Ai are equal, the corresponding matrix is singular. In that case, one may 
use a generalization called confluent Vandermonde matrix [1] [9] , which makes the matrix 
non-singular, while retaining most properties. If Ai = Ai+i = ... = \+k and A.; ^ Ai_i, then 
the (i + fc)th column is given by 
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The confluent Vandermonde matrix looks as 
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Another way to write the (i + k) column is using the derivative, as follows 



(3) 



(4) 



(5) 
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3 The inverse of the Vandermonde matrix 



In applications, a key role is played by the inverse of the Vandermonde and confluent Van- 
dermonde matrices []] [SJ [31 [TOl [II] ?,[T2J[T3]. Both matrices, Vandermonde and confluent 
Vandermonde, can be factored into a lower triangular matrix L' and an upper triangular 
matrix U' where V or C is equal to L'U'. The factorization is unique if no row or column 
interchanges are made and if it is specified that the diagonal elements of V are unity. 
Then, we can write V' 1 = ([/') (L')~ • Denoting ([/') as U, we have found that U is 
an upper triangular matrix whose elements are 

Ui,j = if i > j 

3 1 
Ui.j — J J — otherwise. 

fc=l, k=£i Aj — Afc 

The matrix U can be decomposed as the product of a diagonal matrix D and other up- 
per triangular matrix W. It is very easy to find that 



{ fe=i, fc^ 4 - A fc J (6) 
i+3 



and 

i > j 

W ^ = { n ( A « - X k) otherwise ' ^ 

k—j + l, k^i 

The matrix L — [L') 1 is a lower triangular matrix, whose elements are 

i < j 

Li,j = { 1 i = j . (8) 

Li-ij-i - Lj-ijAi-i i = 2, 3, N; j = 2, 3, i - 1 

Summarizing, the inverse of the Vandermonde matrix can be written as V^ 1 — DWL. 

4 The inverse of the confluent Vandermonde matrix 

We will treat now the case of the confluent Vandermonde matrix. We suppose that just one of 
the values Ai is repeated, and it is repeated m times. We make the usual LU decomposition, 
getting C = L' C U' C1 where L' c is a lower triangular matrix and U' c an upper triangular matrix 
U'. Then, we can write C- 1 = {U' C Y (L' c y . Denoting {U' c )~ as U c , we have found that 
U c is an upper triangular matrix whose elements are 

(U e ) itj =0 i i>3 (9) 

s 

{Uc)ij = tt^tj i = 1,2,3,..., m; j = 1, 2, 3, m (10) 
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{Uc) ^ = ~ji~Ty. S II (Aq - Ag) 

i = 1,2, 3, m; 
j = m + l,m + 2, ...,JV 



(11) 



J 1 

(^)ij= EI n \~T i = m + l,m + 2,...,JV; j = i,...,N (12) 

where it is understood that A m = A TO _i = ... = A2 = Ai, and where the numbers 
A m , A TO _i, A2 appear, they must be substituted by Ai. 

The matrix L c — (L'^ 1 is a lower triangular matrix, whose elements are given by the 
following recurrence relation, 

ifi < j 

(£ck, = { 1 ifi = 3 , (13) 
{L c )i,ji-i,j-i - (L c )i,ji-i,j^i-i ifi = 2, 3, N; j = 2, 3, i — 1 

also here it is understood that A m = \ m -i = ■■■ = A 2 = Ai, and where the numbers 
A TO , A TO _i, A 2 appear, they must be substituted by Ai. 

When more than one value is repeated, the inverse has blocks with the same structure that 
we have already found. 



5 Conclusions 

We have shown a form to determine the inverse of Vandermonde and confluent Vandermonde 
matrices. Although several studies exist for Vandermonde matrices, it is not so for systems 
with repeated eigenvalues, which lead to confluent matrices. Such inverse matrices are of 
importance in several quantum mechanical topics where it is needed to find functions of 
matrices, such as in quantum information processes, where entropies play a key role. 
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